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Critical measurement performance of 
fluid flowmeters requires proper and 
quantified verification data. These data 
should be generated using calibration 
and traceability techniques established 
for these verification purposes. In these 
calibration techniques, the calibration 
facility should be well-characterized and 
its components and performance prop- 
erly traced to pertinent higher stan- 
dards. The use of this calibrator to 
calibrate flowmeters should be appropri- 
ately established and the manner in 
which the calibrated flowmeter is used 
should be specified in accord with the 
conditions of the calibration. 

These three steps: 1) characterizing 
the calibration facility itself, 2) using the 
characterized facility to calibrate a 
flowmeter, and 3) using the calibrated 



flowmeter to make a measurement are 
described and the pertinent equations 
are given for an encoded-stroke, piston 
displacement-type calibrator and a 
pulsed output flowmeter. It is con- 
cluded that, given these equations and 
proper instrumentation of this type of 
calibrator, very high levels of perifor- 
mance can be attained and, in turn, 
these can be used to achieve high fluid 
flow rate measurement accuracy with 
pulsed output flowmeters. 
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Glossarj 


r 




Symbol 


Description 


Dimensions 


Symbol 


Description Dimensions 


Et 


Modulus of elasticity of 


F/V 


Aco.Ac 


Calibrator cross-sectional 
areas at reference and 


V 




the fluid at pertinent 
conditions. 






non-reference conditions, 




Em 


Modulus of elasticity of 


F/V 




respectively. Over-bars, where 






the material of the body 






used, denote averages over 






of the flowmeter at 






piston stroke length. 






pertinent conditions. 




AtfijAi 


Cross-sectional areas of 
tubes connected to 
calibrator piston at 
reference and non- 
reference conditions, 
respectively. 


V 


Et 


Modulus of elasticity of 
the material of the tube 
connected to the 
calibrator piston at 
pertinent conditions. 


F/V 


Co,C 


Circumference of cylindri- 
cal body in unstressed and 
stressed conditions, 
respectively. 


L 


/eo,/e 


Frequency of pulses from 
encoder at reference and 
non-reference conditions, 
respectively. 


pulses/t 



do, d Diameter of rod connected 

to calibrator piston at 
reference and non-reference 
conditions, respectively. 

Do, D Inner diameter of calibrator 
cylinder at reference and 
non-reference conditions, 
respectively. 

Duo, Dm Inside diameter of flowmeter 
at reference and non- 
reference conditions, 
respectively. 

duo, diA Outside diameter of turbine 
propeller at reference and 
non-reference conditions, 
respectively. 

Dca, Dc Inside diameter of the 

cylinder of the calibrator 
at reference and non- 
reference conditions, 
respectively. 

Dm, Dt Outside diameter of the 
tube connected to the 
calibrator piston at 
reference and non- 
reference conditions, 
respectively. 

Ec Modulus of elasticity of the 

material of the calibrator 
cylinder at pertinent 
conditions.. 



F/V 



/mo,/m Frequencies of pulses pulses/t 

from the flowmeter at 
reference and non- 
reference conditions 
respectively. 

Kejh, Ke Encoder constant at pulses/L 

reference and non- 
reference conditions, 
respectively. 

^co, ^c Calibrator constant at pulses/L' 

reference and non- 
reference conditions, 
respectively. 

Kc Calibrator constant in pulses/L' 

pulses per fluid volume 
collected in collection 
vessel. 

^Mo,^M Meter factor at reference pulses/L' 
and non-reference 
conditions, respectively. 
Where over-bars are used 
is meant the average 
value over a Reynolds 
number range. 

Lb),Le Piston stroke length at L 

reference and non- 
reference conditions, 
respectively. 
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Nb>,Ne 


Totalized encoder pulses 
at reference and non- 
reference conditions, 
respectively. 


pulses 


Nm,NM 


Totalized meter pulses at 
reference and non- 
reference conditions, 
respectively. 


pulses 


n 


Unit vector. 


dimensionless 


Pco,Pc 


Pressure of the fluid in 
the cylinder of the 
calibrator at reference 
and non-reference 
conditions, respectively. 


F/L^ 


PcOLLT 


Pressure of fluid in 
collection vessel. 


F/L^ 


Pmo,Pm 


Pressure of the fluid in 


F/V- 



the flowmeter at reference 
and non-reference con- 
ditions, respectively. 

Re Reynolds number. dimensionless 

Ro Roshko number. dimensionless 

5/ Strouhal number. Where dimensionless 

over-bars are used is 
meant the average value 
over a Reynolds number 
range. 

5 Surface area of control 

volume. 

Ti Temperature of meter 

rotor or of rod connected 
to calibrator piston. 

Tco, Tc Temperature of the 

cylinder of the calibrator 
at reference and non- 
reference conditions, 
respectively. 

Teo, Te Temperature of the "C ("F) 

encoder at reference 
and non-reference 
conditions, respectively. 

tco Thickness of the cylinder L 

of the calibrator at 
reference conditions. 

t-m Thickness of the tube L 

connected to the to the 
calibrator piston at 
reference conditions. 



°C ("F) 



ICOLL'T 



Tmo, Tm 



tuo 



u, 



MO 



Fco.Fc 



V 

V 
Fco, Kco 



FcOIX'TO, 
KcOLL'TO 



Fc.Fc 



Fmo> I^mo 



Vm, Vm 



at 



ot 



Time. t 

Temperature of fluid "C (T) 

in collection vessel. 

Temperature of the "C ("F) 

flowmeter at reference 
and non-reference 
conditions, respectively. 

Thickness of the body of L 

the flowmeter at reference 
conditions. 

Average velocity of fluid L/t 

through the flowmeter at 
reference conditions. 

Calibrator piston L^ 

displacement at reference 
and non-reference 
conditions, respectively. 

Fluid velocity vector. 
Control volume. 

Fluid volumetric flowrate 
and volume, respectively, 
in the calibrator cylinder 
at reference condition. 

Fluid volumetric flowrate L7t,L- 
and volume, respectively, 
into collection vessel at 
reference conditions. 

Fluid volumetric flowrate lJ/t,lJ 

and volume, respectively, 

in the calibrator cylinder 

at non-reference 

conditions. 

Fluid volumetric flowrate L^/t,L^ 
and volume, respectively, 
in flowmeter at reference 
conditions. 

Fluid volumetric flowrate lJ/t,lJ 
and volume, respectively, 
through the flowmeter 
at non-reference 
conditions. 

Linear expansion coef- "C"' ("F"') 
of the encoder at 
pertinent conditions. 

Linear expansion coef- "C"' ("F"') 

ficient for the material 
of the calibrator cylinder 
at pertinent conditions. 



yt 

L^ 
V/t,V 



3/*t3 
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Of 




Linear expansion 
coefficient of fluid. 


°C-' fF-^) 


OM 




Linear expansion 
coefficient for the 
material of the body 
of the flowmeter. 


OQ-l (op-lj 


ttr 




Linear expansion 
coefficient for the 
material of rotor or 
turbine wheel. 


"C"' ("F"') 


J8 




Ratio of turbine 
propeller diameter to 
inside diameter of meter. 


dimensionless 


£9 




Hoop strain in (^lindrical 
body. 


dimensionless 


p 




Fluid density. 


MIV 


Pco, 


PC 


Fluid density in calibrator 
at reference and non- 
reference conditions, 
respectively. 


unJ 


PCOLL'TO, 


Density of fluid collected 


wv 


pCOIXT 


in collection vessel at 





Pmh, pu 



Ca 



H), V 



Ma, At 



reference and non- 
reference conditions, 
respectively. 

Density of the fluid in the 
flowmeter at reference 
and non-reference 
conditions, respectively. 

Hoop stress in cylindrical 
body. 

Fluid kinematic viscosity 
at reference and non- 
reference conditions, 
respectively. 

Fluid absolute viscosity 
at reference and non- 
reference conditions, 
respectively. 



M/L^ 

F/V 
V/t 

FtlV 



1. Introduction 

Fluid flowmeter calibration methods use a variety 
of techniques; they include wide ranges of opera- 
tional parameters; they cover wide ranges of preci- 
sion and accuracy, [1-3]. Increasingly, the diverse 
improvements being sought for flowmeters are 



producing corresponding improvements in the 
characteristics of flowmeter calibration systems. 
Of these systems, the piston displacement-type 
calibrator offers advantages such as compactness, 
mobility, efficient change of fluid, and prospects 
for state-of-the-art performance characteristics. 

To characterize the performance of a piston dis- 
placement calibrator, which produces a pulsed 
output signal that is proportional to the volumetric 
flowrate, the objective is, generally, to determine 
the "pulses per volume displacement" ratio (or its 
reciprocal) where the pulse output is assumed to 
come from a source such as a linear encoder. This 
pulse output is also assumed to be proportional to 
the piston displacement. High accuracy calibrator 
performance requires examination of these 
assumptions. 

A number of techniques can be used to deter- 
mine volumetric displacement. It is assumed in 
what follows that the displaced fluid is a liquid, but 
the principles apply to gases or mixtures of gases 
and liquids as well. It is also assumed that both 
temperature and pressure effects on all the compo- 
nents of the piston displacement system, the cylin- 
der, encoder, and fluid should be considered in 
order that high accuracy performance can be 
achieved. 

Depending upon the desired uncertainty level 
for the performance of the calibrator, one or more 
of the pressure and temperature effects on the 
components of the system may be negligible. When 
this is so, it may be permissible to disregard such 
effects to simplify data processing or to reduce the 
size of the controlling software for the system. 
Alternatively, and more preferably, all effects can 
be included in computer software. In this way, the 
terms which are negligible will not influence the 
results when higher levels of uncertainty, i,e,, less 
precise performance can be tolerated or is desir- 
able from benefits vs, costs perspectives. More 
importantly, where high accuracy is required, more 
of the figures available via the software capabilities 
can be accepted as significant. 

The volumetric-type calibrator system, using 
encoded piston displacement as both the flow 
source and as the flow determination scheme, is 
sketched in Fig. 1. The piston motion produces and 
measures a fluid volumetric flowrate that is propor- 
tional to the encoder frequency. The piston in the 
cylinder has a seal that is assumed to seal perfectly 
for all piston velocities. The corresponding fluid 
flowrate through the meter and the meter output 
frequency enable a calibration of the flowmeter. 
Using these elements, three steps are considered. 
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>. Meter pulses 




Piston 
stroke 



■Encoder 
pulses 



Calibrator 



Fig. 1. Sketch of encoded-stroke calibrator in operation. 

The first step is determining the calibrator factor 
which is an "encoder pulses per volume displaced" 
ratio (or its reciprocal) at the defined set of refer- 
ence conditions. These results can be obtained 
experimentally in a number of ways: (1) by physical 
measurement, (2) by the so-called "draw" tech- 
nique, or (3) by using a transfer standard such as a 
single or, preferably, a tandem arrangement of 
calibrated turbine flowmeters. 

The second step is the use of the characterized 
calibrator to calibrate a pulsed-output flowmeter 
such as a turbine meter or other device where the 
meter factor is a "pulses per volume" quantity or 
its reciprocal. This process is done with one or 
more selected fluids and at specified conditions. In 
what follows turbine flowmeters will be considered. 
Conventionally, turbine meter results are produced 
in the form of a meter factor which has units of 
pulses per volume or volume per pulses (refer- 
enced to specified conditions). This meter factor is 
determined over the desired ranges of fluid condi- 
tions and flowrate expressed in terms of a ratio of 
inertial-to-viscous effects such as a Reynolds 
number, or equivalent parameter. 

The third step is the use of the characterized 
turbine flowmeter to calculate a fluid flowrate 
under actual conditions of use. The results can be 
produced with respect to specified reference condi- 
tions or to the actual conditions, depending upon 
the needs of the meter operator. 

The purpose of this paper is to describe these 
three steps and give the pertinent relationships 
that pertain to each procedure. The resulting equa- 
tions are intended to be used in the software pack- 
ages used with these types of calibrators and 
metering units. In this way, it is expected that the 
measurement performance of both the calibrators 
and the metering units can be maximized. 



2. Calibrator Characterization 
2.1 Geometrical Determination at Reference 
Conditions 

To perform the required measurements at refer- 
ence conditions and then calculate the calibrator 
factor in units of pulses per fluid volume displaced, 
the system analyzed is that sketched in Fig. 2. The 
assumption of reference conditions is a conceptual 
situation that is impractical to achieve precisely but 
is done solely for convenience, as will be clear in 
what follows. For the conditions selected to be the 
reference conditions which are denoted by the "0" 
subscript, a specific piston stroke produces a dis- 
placement volume, Vco, and the encoder produces 
the corresponding number of pulses, A^eo. A list of 
symbols is given in the Glossary. The reference 
conditions of calibrator temperature and the pres- 
sure in the calibrator are Tco and i'co, respectively. 
These properties are assumed to be constant and 
steady in the calibrator volume. These conditions 
should be monitored, quantified, and assessed with 
respect to the performance level of the calibrator. 
Furthermore, the fluid inside and outside the cali- 
brator is assumed to have the same temperature as 
the cylinder of the calibrator and no heat is being 
transferred to or from the calibrator. In all that 
follows, the reference conditions of To and Po will 
be assumed to be the same for all components. 




Fig. 2. Encoded-stroke calibrator configuration at reference 
conditions. 



The specified piston stroke, Leo, produces the 
pulse total, A^Eo, where 

NEa=LBaKEo, (1) 

and Keo is the encoder constant in pulses per length 
at the reference temperature condition. To. The 
calibration constant can be written 



Nbj _ Leo Keu Keo 



V —ilSi—tll 

y^ AcoL 



(2) 



COi-'EO 



lea 
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in units of pulses_per volume, where, at the refer- 
ence conditions, Aco is the averaged cylinder cross- 
sectional area over the piston stroke, Leo. This 
calibrator constant, Kco is assumed to be constant 
over the operating range of the calibrator. 

For the calibrator configuration shown in 
Fig. 1, the area Aco is an annular one between the 
cylinder of the calibrator and the rod or tube 
attached to the piston. The precision with which 
^co is determined can be written using root-sum- 
square combinations of component precisions: 



A^Tco , r/AJVEoV , MFcoV]'^ 
Kco [\Nm} \Vco} J 



'ii^hmr- (3) 



where the numerators of the respective terms refer 
to the maximum errors of each of the component 
measurements. Equation (3) indicates that high 
levels of precision in Kco can be attained when 
large pulse sums, ^Veo, and large displaced volumes, 
Fco, are used. Correspondingly, these precision 
levels can be achieved with accurate and sensitive 
linear encoders and accurately measured and large 
cross-sectional areas. 

2.2 Geometrical Determination at Non-Reference 
Conditions 

To perform the required measurements at non- 
reference temperature and pressure conditions, it 
is assumed that these conditions are constant and 
steady as shown in Fig. 3. For a specified piston 
stroke, Le, the corresponding encoder pulse total, 
Ne, is 



Ne=LeKe, 



(4) 



where the encoder constant, Ke, is assumed to 
depend only on temperature according to 



Ke = Keo [1 — oce(Te — Tbo) , 



(5) 



where an is the pertinent linear expansion coeffi- 
cient for the encoder and 7e and Teo are, respec- 
tively, the encoder temperatures at non-reference 
and reference conditions. 

The cross-sectional area change of the calibra- 
tor cylinder that is produced by temperature, Tc, 



and internal fluid pressure, Pc, in excess of the ref- 
erence values is analyzed in Appendix A. In what 
follows, it is assumed that the annular area con- 
tained between the calibrator cylinder and the rod 
or tube connected to the piston is changed with 
temperature in the usual way and the pressure 
effect is considered to enlarge only the cylinder of 
the calibrator. The pressure effect on the tube or 
rod is taken to be negligible. The cross-sectional 
area of the cylinder of the calibrator averaged over 
the stroke length is given, to first order in tempera- 
ture and pressure separately, by 



Ac=Aco[l + 2ac 



iTc-To)][l + 



iPc-Po) Dco 
tcoEf 



],(6) 



where etc is the linear expansion coefficient for the 
material of the calibrator cylinder, Deo and tea are, 
respectively, the inside diameter and wall thickness 
of the calibrator cylinder at reference conditions 
and Ec is the modulus of elasticity of the material 
of the calibrator cylinder. In Eq. (6) the compres- 
sive effects of the pressure in the cylinder on the 
rod attached to the piston are assumed negligible. 
If it is assumed that the linear expansion coeffi- 
cient of the cylinder material is 2 x 10"^ °C"', a sys- 
tematic temperature difference of only 2.5 °C would 
produce a corresponding systematic error of 
±0.01% in the calibrator area. Similarly, if it is as- 
sumed that the modulus of elasticity for the calibra- 
tor material is 2x10 "Pa (2xlO*atm) and if the 
calibrator has a diameter to thickness ratio of 20, a 
systematic pressure error of only 1 MPa (10 atm) 
will produce a systematic error of ±0.01% in the 
calibrator area. Such calibrator area errors will 
propagate through all of the successive relation- 
ships for calibrating and using meters, and these re- 
lationships will be further altered by additional 
temperature and pressure effects. 




i 



Fig. 3. Encodcd-stroke calibrator configuration at non-reference 
conditions. 
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The calibrator constant, Kc, at non-reference 
conditions in units of pulses per volume can be 
written 






(7) 



since Vc is the displacement at non-reference con- 
ditions. Therefore 



Kc=^. 



(8) 



Combining Eq. (8) with Eqs. (2), (5), and (6) pro- 
duces 



Kr = - 



A:co[i-aE(rE-ro)] 



L tcatic J 



or, to first order approximation: 



. (9) 



Kca=Kc 



[l + aniTE-To)] 



[i-2«c(rc-ro][i-fe^] 



.(10) 



For temperatures and pressures higher than refer- 
ence conditions the calibrator constant, Kc, is less 
than the value at reference conditions. In both of 
the above geometrical methods for obtaining cali- 
brator constants, fluid properties are not involved. 
The value for Kco given by Eqs. (2) or (10) 
should be put into the software that operates the 
displacement calibrator, together with pertinent 
material constants and component dimensions. 
Then the computation for the calibrator constant 
Kc at non-reference conditions can be done using 
appropriate measurements for encoder tempera- 
ture, in addition to the pressure and temperature 
in the calibrator, as shown in Eq. (9). These values 
of Kca or Kc will be used for the accurate calibra- 
tions of flowmeters. 

2.3 Determination of Calibrator Constant by the 
Draw Technique at Reference Conditions 

The configuration for the draw technique is 
shown in Fig. 4, where the piston stroke displaces 
an amount of fluid through a valve into the collec- 
tion tank; the corresponding encoder pulses total 
Neo. All temperatures and pressures are the refer- 
ence values. In all that follows, it is assumed that 



the seal of the piston in the calibrator seals per- 
fectly and is sufficiently pliant to continue to seal 
perfectly when changes in temperature and inter- 
nal pressure occur and change the diameter of the 
cylinder of the calibrator. Applying conservation of 
mass principles to the constant volume shown in 
figure 4, we obtain: 



I//-"--! 



pvn dS =0, 



(11) 



where p is the fluid density, V is the control volume 
surrounded by the control surface, S,v is the fluid 
velocity vector, and n is the unit vector normal to S 
with positive direction pointing to the interior of V. 
When 



ll/'^»^=^' 



(12) 



there is no change of the mass within the control 
volume in time. This means the effects of fluid 
friction or of heat transfer are negligible or, taken 
in total, do not alter in time the mass contained 
within the control volume. Under such conditions, 
the inlet and outlet mass fluxes through the control 
surface, S are the same, that is 



pea Kco — PCOLUTO KcOLL'TOj 



(13) 



where pco and poollto are the fluid densities, 
respectively, in the calibrator and collection tank at 
the same reference conditions and therefore are 
equal. The quantities Fco and Kcoll-to are the volu- 
metric flowrates out of the calibrator and into the 
collection tank, respectively. 



Control volume, V 
surrounded by 
control surface, S 

Valve 




Vcoll'to 

TcoD'to 



Hg. 4. Experimental arrangement for the draw technique at 
reference conditions. 



Therefore, since the times of displacement and 
collection are assumed to be the same. 



Kco=K( 



COLUTO > 



(14) 
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and the calibrator constant, in units of pulses per 
displacement volume can be written 



where 



Kai=Nm/VcouLTo, (15) 

Vcoix"n= AcoLeo- (16) 

Therefore, using Eq. (13)-(16), we obtain 

Ka,=Km/Aco, (17) 

in agreement with Eq. (2). Accordingly, results 
obtained via this procedure should duplicate those 
obtained via the geometrical measurement 
techniques. 

It is generally assumed that the calibrator con- 
stant Kco is independent of the piston speed or the 
fluid properties. If this assumption is not valid then 
it is necessary to characterize the dependence. This 
should be done in a manner analogous to the pro- 
cedures used to characterize turbine meters as 
given below. 

It should be noted that, for the control volume 
shown in Fig. 4, if 



iJ/^^-'O' 



(18) 



then the mass flux displaced in the cylinder of the 
calibrator is not equal to that delivered to the col- 
lection tank. Furthermore, if the fluid density in 
the control volume decreases because of a temper- 
ature rise then the mass flux to the collection tank 
would be more than that displaced in the calibrator 
cylinder. As will be seen below, this produces a 
turbine-type meter calibration result in pulses per 
volume units that is larger than it should be. Con- 
versely, if the fluid density in the control volume 
increases, then the mass fliix delivered to the col- 
lection vessel is less than that displaced in the cali- 
brator cylinder. This produces a low calibration 
result for a turbine meter constant in pulses per 
volume units. 

2.4 Determination of Calibrator Constant by the 
Draw Technique at Non-Reference Conditions 

The configuration for the draw technique in non- 
reference conditions is shown in Fig. 5. As before, 
the piston stroke displaces a volume of fluid in the 
cylinder and the corresponding encoder pulses 
total Ne. The temperature and pressure of the fluid 
in the cylinder are Tc and Pc, respectively. Apply- 



fc.Np 



Piston 

stroke 



D 



Vc 




Control vclume.V 
sun'oundsd by 
conJnol surface, S 

Valve 



Vooll't 

Tooirt 
Pcoirt 



Fig. 5. Experimental arrangement for the draw technique at 
non-reference conditions. 



ing conservation of mass principles to the control 
volume in Fig. 5, for the condition where there is no 
change with time of the mass within the control 
volume, Eq. (11) and (12) indicate that 



pc Vc= Pcoll'tVcolitt , 



(19) 



or when the displacement and collected volumes 
are simultaneous 



pc Vc= PcoivtVcoll't > 



(20) 



where pc and pcoll't are, respectively, the density of 
the fluid displaced within the cylinder of the cali- 
brator and flowing from the valve and into the col- 
lection tank. The volume Kcoll-t can be determined 
using volumetric or gravimetric techniques, or both. 
It is noted that constant conditions of temperature 
and pressure need to be maintained in the calibra- 
tor and in the collection vessel for accurate results. 
With the volume Fcollt determined, Vc can be 
computed via 



Vc= {pCOUJil pc) VcOLVT, 



(21) 



where, according to the definitions of fluid thermal 
expansion effects and compressibility, to first order 
in temperature and pressure, separately: 



PCQLLT ^ 

Pc 

[1— 3aF(TcOLL'T 



-Tc)-\[l + 



(/*COLL'T— -Pc)] 



(22) 



where of is the linear expansion coefficient for the 
fluid, the pressures Pcollt and Pc are, respectively, 
those for the fluid passing from the valve into the 
collection tank and in the cylinder of the calibrator, 
and Ef is the modulus of elasticity of the fluid. It is 
noted that both m and Ev are dependent upon 
temperature and pressure for the specific fluid, but 



516 



Volume 97, Number 5, September-October 1992 

Journal of Research of the National Institute of Standards and Technology 



the values used here and in what follows are 
assumed to be averages taken over the appropriate 
ranges of temperature and pressure. It is also 
noted here that the reciprocal of the modulus of 
elasticity of the fluid is also the compressibility of 
the fluid. 

In Eq. (22), it is noted that temperature and 
pressure effects have opposite signs in producing 
fluid density changes. However, when the tempera- 
ture differences between the collection vessel and 
the calibrator are large and where the pressure in 
the calibrator is much larger than that in the 
collection vessel, the effects on the fluid can be 
signiflcant. 

If it is assumed that the linear expansion coeffi- 
cient of a hydrocarbon liquid is 3 x 10"'* °C~', a sys- 
tematic temperature error of only 1 °C will produce 
a systematic error of 0.1% due to temperature in 
the determination of the calibrator volume by this 
draw method. Similarly, if it is assumed that the 
fluid's modulus of elasticity is 2 x 10' Pa (2 x 10'' 
atm), a systematic pressure error of only 2 MPa 
(20 atm) will produce a systematic error of ±0.1% 
in the calibrator volume. As mentioned above, such 
errors will propagate through all of the relation- 
ships for calibrating and using meters and these 
errors will be further altered by additional temper- 
ature and pressure corrections. 
Furthermore, using Eq. (7) we can write 



Kc = 



Ne 



/ PC0LL'T ^ 



(23) 



Vo 



OLL'T 



then, from Eq. (22) 



^E 



Kc = rr^^- [l + 3ar (Tcolvt - Tc)] 

KCOIX'T 

^|-^_(PcoigpPc)j 



(24) 



This result is then related to reference conditions 
via Eq. (10) to give 

^CO" r (p _p ^-l 

Ne [1 + 3aF(rcoLL'T- Tc)][ 1 - """^^ p + "eCTe- T„)] 



Kcou.T[l-2ac(rc-rn)] 1 -— r . 

L ICOiic J 



(25) 



It is noted that the cross-sectional area, Ac can be 
determined via Eqs. (8) and (23) to be 



^'^-ateI PC r^"^""- 



(26) 



Using Eqs. (5), (6) and (22), this can be written in 
terms of directly measured quantities and related 
to reference conditions via 



/f Eo P'cou.'T [1 - MTe- To)] [1 - 2ac (Tc- r„)][l -^^y^^] 



^E[l+3MWT-rc)][l-^^^!g=^] 



If Kc is defined as 



K.c—NElV<xi\.irt , 



(27) 



(28) 



then 



A CO ~ '^C X 



{Pc-Pa)Do3-\ 



{\+2ac(Jc-m[\+——^ J 



[1 -3aF{rcoLLT- rc)][l +££2I£IZ£^ j|;i _ a^iJn-T,)] 

or, to first order 



(29) 



{Pc-Po)Da,-\ 



[l-2ac{Tc-T„)] 1 — 

L tcotc J 



[1 +3«F{7'coa'T- rc)][i -(^£2!|r:^ ][i+aE(rE-r„)] 



.(30) 



When temperatures exceed those of the reference 
conditions and pressures are those for the refer- 
ence conditions: 

K6<Kco. (31) 



Thus, the observation made above is repeated 
here, namely, that when pressure effects can be 
neglected and when temperatures are above the 
reference conditions, the calibrator delivers, for 
the same encoder output pulses, more fluid volume 
than would occur under reference conditions. 



517 



Volume 97, Number 5, September-October 1992 

Journal of Research of the National Institute of Standards and Technology 



It is also noted that Eqs. (23)-(27) can be 
combined to give: 



Kc=Kco~ 



[l-asiTE-n)] 



[l+2ac(Tc-To][l+^^^=^^] 



,(32) 



which is the same as Eq. (9). 

It is noted that, in these two draw procedures, 
the decision to collect a sufficiently large number 
of encoder pulses should precede the operation. 
This number should be selected according to the 
desired precision for the calibrator constant, see 
Eq. (3). In the reference conditions 

Nb>=Kb,Leo, (33) 

and in the non-reference conditions 

Ne=KeLe, (34) 

but, to first order in temperature, 

Ke =Kw[1- aE (Te - To)} , (35) 

and 



LE=LEo[l + aE(7'E-ro)]. 



(36) 



Combining Eqs. (33)-(36) therefore yields: 

NE=KELE=KB>LBi=NEo. (37) 

This indicates that the precision criterion 
specified for the calibrator should be achieved via 
the number of pulses selected and the choice is 
not dependent upon whether reference or non- 
reference conditions prevail. 

As stated above following equation (17), the cali- 
brator constants determined via the draw proce- 
dure are generally assumed to be independent of 
piston speed and fluid properties. Where this is not 
valid, efforts should be made to achieve this 
assumption, i.e., improving the piston seals or the 
calibrator should be characterized using techniques 
analogous to those for turbine meters as will be 
described below. 

It should also be noted that in the above 
described draw procedures, the valve and pulse 
counting techniques must not introduce spurious 
effects. Spurious counting effects may result from 
the fluid dynamics in the valve as the flow is started 
and stopped in conjunction with starting and stop- 



ping the pulse count. If such an effect is present 
then it should be eliminated or proper account 
made for it so that the appropriate fluid volume is 
associated with the pulse total. This volume can be 
determined using the valve compensation tech- 
niques that are conventionally applied to diverter 
systems, see [4]. 

To use the now-characterized calibrator, it is 
required that appropriate instrumentation be prop- 
erly installed both to assure that Eq. (12) is satis- 
fied and to measure the quantities involved in Eqs. 
(9) or (10) or (32) for the calibrator and equations 
(28) or (29) or (30) for dispensing precise volumes 
of fluid at specified conditions of temperature and 
pressure. Appropriate values are needed for the 
material constants— the thermal expansion coeffi- 
cients aE, ap and ac, the modulus of elasticity of 
the material of the calibrator cylinder, Ec, and that 
of the fluid, Er, and the pertinent dimensions of 
the cylinder. Once the appropriate value for the 
constant Kco has been installed in the calibrator 
software or the working procedures for the cali- 
brator, the next step is to use the calibrator to 
calibrate a flowmeter. 



3. Calibrator Use in Calibrating a 

Turbine Type Flowmeter 
3.1 Reference Conditions 

To calibrate a turbine-type flowmeter at refer- 
ence temperature and pressure conditions using 
the calibrator characterized as described above, the 
arrangement is sketched in Fig. 6. The temperature 
and pressure are the reference conditions denoted 
by To and A. As stated above, ro=rMo=7'co and 
Po=Pmo=Pco. Applying again the conservation of 
mass principles, . and assuming that there is no 
change of mass within the control volume with 
time, Eqs. (11) and (12) indicate that 



pco f^co — pMo Vmo , 



(38) 



where pco and pMo are the fluid densities in the 
calibrator and meter, respectively, at the reference 
conditions. The quantities Kco and Vm are the volu- 
metric flowrates, respectively, in the calibrator and 
through the meter at the reference conditions. 

Since reference conditions in cylinder and meter 
are assumed the same, the densities pco and pMo 
are equal and therefore: 



Ka) = KMo. 



(39) 
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Fig. 6. Arrangement for using characterized calibrator to cali- 
brate a turbine-type flowmeter at reference conditions. 

Since the pertinent time intervals are assumed the 
same: 



^EO NeB' /eo 

Kco=-:z~-~z. "=~i~t 

A CO A CO Leo Kco 
where /eo is the encoder frequency, and 



(40) 



A'^MO /mo 
^M0= =T— , (41) 

f^MO I^MO 

where /mo is the meter frequency. 

In Eqs. (40) and (41), the pulse counts A^eo and 
Nmo are to correspond to the respective volumes Vco 
and Fmo which are assumed to be the same, or 
appropriate corrections are made to compensate 
for any differences between these volumes. 
Combining Eqs. (40) and (41) gives 

/mo -Keo _Nmo Keo _A^mo _/mo 



/eo Aco ^eo Aco ^e» 



/e. 



(42) 



where ^mo has units of pulses per volume at 
the reference conditions. The quantity Keo is ob- 
tained from manufacturer's specifications, testing 
or Eq. (1). 

Conventional non-dimensionalization procedures 
can be applied to the flowmeter characteristics to 
produce a set of parameters which interrelate the 
significant inertial, viscous, and oscillatory effects 
that constitute the performance of the meter in the 
calibration conditions, see [6,7]. In this way, the 
performance of the device can be predicted for 
other fluid and flow conditions where this set of 
parameters are the pertinent ones to describe the 
meter's performance. Of course, when other, differ- 
ent effects such as fluid compressibility or gravita- 
tional influences become significant, it should be 
expected that the initial parameterization needs to 
be modified to include such effects to obtain satis- 
factory description of meter performance, see [8,9]. 



By normalizing the meter factor, Kuo and the 
fluid flowrate using the^ meter diameter, Dmo and 
the mean flow velocity FmoMmo= Umo and the fluid 
kinematic viscosity vo, we can obtain, for example, 
the Strouhal number 

St = — — = Ci KmoD mo « KmoDmb , (43) 

[/mo 

where Ci =it/4. The Strouhal number is the ratio of 
characteristic meter frequency effects such as 
propeller rotation rate to fluid momentum effects. 
By effects here is meant either forces or energies. 
As such, the Strouhal number is a dimensionless 
meter factor. 

Conventionally, it is fluid mechanical practice to 
formulate the Reynolds number as the ratio of 
fluid inertial to viscous effects. 



Re = 



D moUmo 

Vo 



(44) 



This is generally used to describe the domain of the 
meter calibration for which the corresponding 
Strouhal numbers specify the range of meter 
response. A Strouhal-Reynolds characterization of 
a pulse-producing flowmeter conforms to conven- 
tional fluid mechanical procedures and it is 
analogous to orifice metering practice where 
discharge coefficient (a ratio of fluid inertia to 
differential pressure effects) is described function- 
ally or graphically versus Reynolds number. How- 
ever, as noted in Eq. (44) the Reynolds number 
requires that the fluid velocity be known. Since this 
is the purpose for using the meter, an iteration 
technique is required to calculate the flowrate. To 
avoid such an iteration, turbine meter manufactur- 
ers have designed their products to have high levels 
of linearity over wide flowrate ranges. As a result 
of this, it has been conventional turbine meter 
practice to characterize performance via K-factot 
vs. frequency to kinematic viscosity ratio— the 
so-called Universal Viscosify Curve (UVC). In 
accord with the principles of dimensional simili- 
tude, the meter frequency to fluid viscous effects 
can be formulated using, as characteristic length 
scale, the meter diameter, £> mo as in Eq. (44). This 
formulation can also be achieved via the product of 
Strouhal and Reynolds numbers; this product has 
been recently referred to, see [9], as the Roshko 
number. 



Ro = 



fuoDuo 



Vo 



(45) 
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This Roshko number is the dimensionless 
version of the frequency-to-kinematic viscosity 
ratio used for the UVC. It is expected that the 
dimensionless version of the UVC— the Strouhal- 
Roshko characterization— should produce superior 
prediction of meter performance compared to the 
UVC by virtue of its more complete, i.e., dimension- 
less description of the inertial, viscous, and oscilla- 
tory effects that occur in the meter operation. 

The question of which set of parameters best 
describes a turbine meter's performance should be 
determined using appropriate data sets. The 
specific design considerations of blade size and 
shape, internal meter geometry, bearing features, 
etc., can be expected to play significant roles in the 
selection of non-dimensional parameters. Where it 
may happen that the Reynolds number produces a 
better description of turbine meter performance, 
i.e., the data collapses better onto a single curve 
than achieved using the Roshko number, an itera- 
tive sequence of computations may be necessary to 
produce an accurate determination of flowrate. If 
the meter is very linear, such an iteration may not 
be required, laical flowmeter characterization 
results are sketched in Fig. 7. Because it is not 
known whether Reynolds or Roshko number is the 
better parameter with which to characterize the meter 
performance, both shall be included in what follows. 



St 




Re or Ro 



Fig. 7. Normalized calibration results for a turbine-type flow- 
meter. Points denote averaged results; bars denote standard 
deviations obtained at each flowrate. 



In those instances where meter performance 
might deviate from the curve shown in Fig. 7, the 
interpretation would be that additional factors in 
the deviant conditions have become significant 
whereas these factors were insignificant in the cali- 
bration conditions. Examples might be fluid 
frictional effects in the turbine bearings produced 
by extreme viscosity variations from those prevail 



ing in the calibration conditions, or liquid cavita- 
tion effects, etc. 

The curve shown in Fig. 7 is interpreted as the 
functional relationship between the Strouhal and 
Reynolds or Roshko numbers which are assumed 
to be the salient parameters describing the perfor- 
mance of this flowmeter over these calibration con- 
ditions. Inherent in this interpretation is the 
assumption that in any subsequent use of this func- 
tional relationship, the geometries of this meter, 
for example, the bearings or the propeller diameter 
are not changed relative to the selected character- 
istic length of the meter, i.e., the internal diameter, 
D. If such changes do occur it can be expected that 
the curve shown in Fig. 7 can change. For example, 
if a smaller propeller were installed, this curve may 
retain its shape but lie below that shown in Fig. 7. 
A non-dimensional parameter which could take 
into account different propeller diameters is 
P=dlD, where d is the propeller diameter. 
Accordingly, the curve shown in Fig. 7 would per- 
tain to the specific /3 for which the calibration was 
done. If other, smaller propellers were also cali- 
brated, these results could be plotted in Fig. 7 and 
parameterized with the smaller value of p. When 
the propeller and meter-body materials are the 
same and where pressure effects can be neglected, 
the p ratio will remain constant when the tempera- 
ture changes. When the propeller and the meter- 
body materials are different, temperature changes 
can produce different ^ ratios. These different 
ratios can be computed using pertinent relation- 
ships. The computed results should then be used 
with calibration data taken for different ^ ratios to 
predict the meter performance at the different 
temperature conditions, see Appendix B. 

Figure 7 is different from conventional turbine 
meter performance plots in which the meter's 
A^-factor is plotted versus the ratio of frequency-to- 
kinematic viscosity. These conventional quantities 
are different from the Strouhal and Reynolds or 
Roshko number parameters by constant factors 
and by factors of the meter diameter raised to dif- 
ferent exponents. These powers of the diameter 
should change only slightly with small changes in 
temperature and internal pressure. However, when 
conditions vary widely, the dimensionless formula- 
tions should be used and are expected to produce 
improved meter performance. Plots in the format 
of Fig. 7 should then, with the exceptions of devi- 
ant phenomena becoming influential, apply to a 
wide range of specific, dimensional fluid property 
and flow conditions and produce accurate predic- 
tions of turbine meter performance. 
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Where conventional turbine meter practice is 
used and plots are produced for K-factor versus the 
ratio of frequency-to-kinematic viscosity, improved 
meter performance can be expected when these 
calibration results are corrected to specified refer- 
ence conditions. Accordingly, the corrected results 
should Incorporate temperature and internal pres- 
sure corrections for the ratios of meter diameters 
raised to the relevant exponent to predict meter 
performance for other fluid and flow conditions. 
For meters having good linearity characteristics, 
i.e., constancy of the meter factor over specified 
flowrate ranges, the more important of these two 
corrections is that for the ^-factor, i.e.. 



Kmo—Km {Dm/Di^) 



(46) 



This stipulates that the meter frequency-to-fluid 
inertial effects ratio, i.e., Strouhal number be the 
same in the actual conditions as in the reference 
conditions. To first order in temperature and pres- 
sure, separately, the diametral ratio can be written, 

(#^r = [l+3aM(rM-7'Mo)] 

\iJMO' 

r^ 3(Pm~^Mo) ^MO |_ ^M0 fA~^ 



where Dm and Dmo are, respectively, the meter 
diameters at the non-reference conditions, Tm and 
Pm and reference conditions, Tm and Pm- The 
quantities um and Em are, respectively, the linear 
expansion coefficient and the modulus of elasticity 
of the meter body material for these conditions, 
and tMo is the thickness of the meter body at refer- 
ence conditions. The more significant of the two 
correction factors is usually that for temperature 
deviations from reference conditions. When the 
pressure correction can be neglected. 



Kmo — Km [1 + 3aM (^m — Tmo)] , 



(48) 



or, to first order: 



Km =Kmo [1 - 3aM(rM - Tmo)] . (49) 

This relationship duplicates that given in [10]. 

3.2 Non-Reference Conditions 

To calibrate a turbine-type flowmeter using the 
calibrator characterized as described above, the 
arrangement is sketched in Fig. 8. The temperature 



and pressure in the calibrator are the non- 
reference conditions denoted by Tc and Pc. Apply- 
ing again the conservation of mass principles, and 
assuming that there is no change of mass within 
the control volume with time, Eq. (11) and (12) 
indicate that 



PcVc= PmVm ' 



(50) 



Since the fluid conditions in the cylinder and in the 
meter can be different, the ratio of the densities is 



^=[l-3a^iTc-TM)][l+^^^], (51) 



and from Eq. (32) and in a manner analogous to 
Eq. (40) 

J, Ke /e 



Ac Vc 



K 



[1-aEiTE-To)] 



CO 

[1 +2ac (Tc- To] [l +^^£^^^] 



(52) 



Then, analogously to Eq. (42) 



Am-t-Ac — -Tr-Aco 

/E Pc Ne 



[1 - aE(rE-ro)][i+3aF(rc-rM)] 



(53) 



where ^m has units of pulses per volume at the 
specific, non-reference meter conditions. 



fp ,Np 
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Fig. 8. Arrangement for using characterized calibrator to cali- 
brate a turbine-type flowmeter at non-reference conditions. 
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As stated above after Eq. (41), the pulse totals 
from the meter and the encoder have to corre- 
spond to the same displaced volume or time inter- 
val. If compensations are needed to achieve this 
correspondence, these should be done and results 
used in Eq. (53). 

By normalizing this meter factor, ^m and the 
fluid flowrate using the non-reference meter 
diameter. Dm, and fluid kinematic viscosity, v, we 
obtain Strouhal and Reynolds or Roshko numbers 



St —/mDm/Um—CiKmOm , 



where Ci='k/4, and 



Re =DmUm/v 



or 



Ro = 



/mDm 



(54) 



(55) 



(56) 



Using these dimensionless parameters, the per- 
formance for the meter can be plotted; results 
should be as shown in Fig. 7. As noted above in 
Eqs. (43)-(45), the characteristics of meter factor 
and diameter at reference and non-reference 
conditions are interrelated and the dependence of 
the fluid's kinematic viscosity can be written func- 
tionally as 



v = vo[T,To,P,Po]. 



(5.7) 



With the performance curve given in Fig. 7 and 
the relationships given in Eqs. (43)-(45), or (54)- 
(57), one is now ready to use the flowmeter to 
make a flowrate measurement. 



4. Using a Turbine-T^pe Flowmeter 

To Make a Measurement 
4.1 Reference Conditions 

Given that the meter performance characteristics 
are as shown in Fig. 7, or less preferably but more 
conventionally as in Fig. 9, one can quantify the 
meter linearity over a specified flowrate range. The 
meter linearity is conventionally the average of the 
maximum and minimum values of the meter factor 
(Strouhal number) over this range; normalized by 
the average meter factor, see [3 and 10]. For the 
meter performance shown Jn Fig. 7, the mean value 
of the Strouhal Number, St gives the mean meter 
factor, via 



^MO— " 



St 



CiD 



(58) 



IJ^MO 



where Ci='ir/4. The flowrate is determined via 

.3 



• _ /mo C\Di/iafuo 

J'mo— ^ — — zz • 

Kuo St 



(59) 



This result would pertain to any flowrate over the 
range specified for the meter's linearity. 



K. 




f/v 



Fig. 9. Conventional performance plot for a turbine-type 
flowmeter. 



If it is desired to improve the accuracy of this 
flowrate determination, this could be done by 
using, for example, the curve shown in Fig. 7 or 
close approximations to it. With advances in 
today's computer technology, this type of process 
can be readily installed in the secondary devices 
used with flowmeters. For the specific frequency 
from the meter, the Roshko number can be calcu- 
lated directly and then used to determine the 
corresponding Strouhal number. For the case 
where the meter is characterized using Strouhal 
and Reynolds parameters, the process to deter- 
mine an accurate flowrate should be iterative. This 
iteration process should begin using a mean value 
of meter factor, such as given in Eq. (58), this value 
of meter factor enables a computation of the 
flowrate via Eq. (59). Using this flowrate, the 
Reynolds number can be computed and then used 
to get the corresponding Strouhal number from the 
calibration curve and a refined value of flowrate. 
This process should be repeated until satisfactorily 
small changes are found in successive results. In 
this way, the accuracy level for the flowrate deter- 
mination can be increased over the level associated 
with the meter's linearity. This can be done to the 
precision level associated with a specific flowrate as 
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quantified in the calibration process. The resulting 
enhanced meter performance could amount to 
significant improvements in measurement accu- 
racy. In the following, a Strouhal-Reynolds charac- 
terization of meter performance will be used since 
it is more conventional in fluid mechanics and 
since it may require the iteration procedure, 
described above, to be used. 

4.2 Non-Reference Conditions 

To describe meter performance in non-reference 
conditions, the meter characteristics shown in 
Fig. 7 will be used. The reason for this is that the 
complete, non-dimensional assessment of meter 
frequency effects and fluid inertial and viscous 
effects are not complete in Fig. 9. For this reason, 
the non-reference meter diameter. Dm, and 
frequency, /m, and the fluid's kinematic viscosity, v, 
should be used, in compatible units, to produce the 
Reynolds number for the non-reference conditions. 
This Reynolds number produces, using Fig. 7, the 
corresponding Strouhal number which with the 
meter diameter. Dm, gives the meter factor, ^m. 
The flowrate measurement is then obtained using 



• _/m _Ci/m-Dm 

•^M St 



(60) 



This flowrate is correctly converted to reference 
conditions by specifying that the Strouhal number 
is the same for the non-reference and reference 
conditions, specifled by this Reynolds number, i.e., 



- „ fuDu fuoDm 
ii/=Ci — ; — = Ci — z . 



Km 



Kmo 



(61) 



Therefore 



• • /mo/Dmo\' .,„ 

VMo = VM-r{j^) . (62) 

The reference to non-reference frequency ratio is 
obtained by specifying that Reynolds number 
similarity exists for these two conditions, i.e.. 



Re = 



MO 



/mDm _ 

(5/) wm (5r) vm ' 



fmD 



(63) 



It is noted that this is equivalent to stipulating that 
Roshko number similarity exists for these two 
conditions. From Eq. (63) 



/mo_i'mo /Dm\^ 
/m Um V-Dmo/ 



(64) 



and 



[l+3aF(r„-r„„)][l-^^i!^^] 



) /Hmov 



(65) 



It is apparent that, to obtain high accuracy 
flowrate measurements using the procedures 
described above, appropriately high accuracy 
measurements are required for the component 
measurement systems and for the pertinent mate- 
rial properties as shown in Eq. (65). In turn, it 
appears feasible that once systematic uncertainties 
are satisfactorily removed from calibration 
facilities, the measurement processes in calibration 
laboratories, and the measurement systems making 
on-line measurements will be commensurately 
improved and the uncertainty levels for these 
measurements can be predicted using such Eqs. as 
(2), (10), (25), and (65). 

5. Discussion 

The above-derived results can, for the sake of 
brevity, be assessed by considering the effects of 
temperature and pressure on the respective 
factors — the calibrator and flowmeter constants 
and the flowrate measurement at reference condi- 
tions. To do this Table 1 shows, for specific mate- 
rial conditions and geometrical sizes, the variations 
associated with temperature and pressiire effects 
separately and then summed together. 

The results shown in the first row of Table 1 
present the variations in calibrator constant Kca 
determined using geometrical measurement meth- 
ods for variations of ± 1 °C in temperature and 
±1x10* Pa (1 atm) in pressure. The material 
properties and geometrical assumptions are given 
under the headings of the respective columns. The 
worst-case combination of temperature and 
pressure variations taken separately are given in 
columns 7 and 12, respectively. These results are 
obtained by adding the absolute values of the com- 
ponent contributions. The total worst-case combi- 
nation for temperature and pressure variations 
taken together is given in the column at the right 
side of the table. Accordingly, the total tempera- 
ture effect on Kca in the measurement method is 
five times larger than the pressure effect and the 
total of these gives an imprecision of ± 0.006%. 
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Table 1. Temperature and pressure effects 


















Factor 
affected 


Total 
diff 

CO 


Temperature effects 


Total 

temperature 

combined 

worst case 


Press 
dice 
(atm) 


Pressure effects 


Total 

pressure 

combined 

worst case 


Total 
worst case 


Encoder 


Calibrator Fluid 


Meter 


Calibrator 


Fluid 


Meter 






factor 


cross-section expansion 
') (oc=2xlO-'°C-') (of=3x10-'°C- 


expansion 
') (o„=2xlO-">C-') 






cross-section 
(C/j=20; 
£c=2xlO"Pa) 


expansion 
(£F=2xlO'Pa) 


expansion 

(C/J=15; 

£M=2xlO"Pa) 






if CO : CEOm" 


il 


±0.001% 


±0.004% 




±0.005% 


±1 


±0.001% 






±0.001% 


±0.006% 


i^CO : DRAW 


+ 1 


1:0.001% 


±0.004% ±0.090% 




±0.095% 


+ 1 


±0.001% 


±0.005% 




±0.006% 


±0.101% 


Kmo' 


il 


±0.001% 


±0.004% ±0.090% 




±0.095% 


±1 


±0.001% 


±0.005% 




±0.006% 


±0.101% 


Sf 


+ 1 


±0.001% 


±0.004% ±0.090% 


±0.006% 


±0.101% 


±1 


±0.001% 


±0.005% 


±0.001% 


±0.007% 


±0.108% 


K»' 


^ 1 


±0.002% 


±0.008% ±0.180% 




±0.190% 


±1 


±0.002% 


±0.010% 




±0.012% 


±0.202% 


Si' 


±1 


±0.002% 


±0.008% ±0.180% 


±0.006% 


±0.196% 


± 1 


±0.002% 


±0.010% 


±0.001% 


±0.013% 


±0.209% 


VmJ' 


+ 1 


±0.001% 


±0.004% ±0.090% 




±0.090% 


±1 


±0.001% 


±0.005% 




±0.006% 


±0.101% 


Km" 


£1 


±0.002% 


±0.008% ±0.180% 




±0.190% 


±1 


±0.002% 


±0.010% 




±0.012% 


±0.202% 


Kmo' 


+ 1 


±0.002% 


±0.008% ±0.270% 


±0.012% 


±0.292% 


+ 1 


±0.002% 


±0.015% 


±0.002% 


±0.019% 


±0.311% 



* Determined via Eq. (10) at non-reference conditions. 

** Determined via Eq. (29) at non-reference conditions with no additional uncertainty 

for the encoder or the collection volume. 
' Determined using .ATci): draw via Eq. (42) with no additional uncertainty resuhing 

from encoder or meter outputs. 
"■ Determined using Kco -. draw via Eq. (43) in reference conditions. 
° Determined using ^cii : draw via Eq. (53). 



'Determined using Kca-. draw via Eq. (54). 

^ Determined using ^mi> via Eq. (59) with no additional uncertainty resulting 

from the meter output. 
'■ Determined using Ku via Eq. (60) with no additional uncertainty resulting 

from meter output. 
' Determined using Eq. (65) with not additional uncertainty for meter output, 

reference meter dimension, and absolute viscosity. 
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The results for ^co determined using the draw 
technique are shown in the second row. These give 
a total temperature variation of ±0.095% owing 
mainly to fluid expansion effects. The pressure 
variation is noted to be a factor of about 16 less 
than this level and the total imprecision totals 
±0.101% which is about a factor of 20 larger than 
the level achieved using the geometrical measure- 
ment technique. 

The results for the meter factor, Kmo, are given 
in the third row of the table. Here, at reference 
conditions, no additional uncertainties are shown 
over those for ^co in accord with Eq. (42). 

The fourth row indicates that the total tempera- 
ture variation for St is increased over that for Kuo 
by the amount allocated for meter expansion. This 
produces the ±0.101% uncertainty total for 
temperature which when added to the increased 
pressure total gives ±0.108%. This shows that 
while the Strouhal number is the preferred non- 
dimensional parameter to characterize the fre- 
quency effects of a meter, it has more factors than 
the conventional meter factor, ^mo and, therefore, 
it can have increased uncertainty. However, 
because the Strouhal number is a dimensionless 
ratio of frequency to inertial effects in the meter it 
should be successful in producing more satisfactory 
metering results for widely ranging conditions than 
can be done using a dimensional quantity such as 

The fifth and sixth rows express the uncertainties 
for Km and the corresponding Strouhal number in 
a manner analogous to that used to obtain values in 
rows three and four. However, as shown in Eq. (53) 
the uncertainties attributed to Kco are increased 
appropriately due to encoder, calibrator, and fluid 
effects; and, in the case for St, meter effects. It is 
noted that the totaled uncertainties have essen- 
tially doubled in comparison with those for Km. 

While it is recognized that the conditions of 
± 1 °C and ± 1 X 10^ Pa (1 atm) pressure variation 
can be termed large, it should be recognized that no 
uncertainty has been allocated in Table 1 for the 
uncertainty with which the material properties and 
geometrical dimensions are known or have been 
determined. As well, it should be recognized that 
these assumed conditions of ± 1 °C and ± 1 X 10^ 
Pa (1 atm) can in reality understate the actual vari- 
ations that may exist in real situations where 
flowrate measurements can be attempted in harsh, 
hostile environments. 



It should be emphasized that the above discus- 
sion deals only with the uncertainties associated 
with temperature and pressure imprecision and 
does not take into account contributions from 
flowmeter imprecision such as the variability of 
detecting meter pulses. Additionally, systematic 
errors which can greatly exceed the levels of impre- 
cision are not included in any of the above. 

A summary of the equations derived for charac- 
terizing piston-type, encoded-stroke calibrators, 
their calibration procedures, and the subsequent 
use of calibrated turbine meters is given in 
Appendix C. These equations are found to be quite 
simple when reference conditions prevail. How- 
ever, reference conditions are practically fictitious 
and therefore the performance of all of these 
devices should be considered as occurring in non- 
reference conditions. Given the capabilities avail- 
able in today's computers, it is felt that the 
complete equations for these types of devices 
should be installed in the controlling and process- 
ing software so that when temperature and or pres- 
sure conditions become significant in these 
processes the results are accurate. Using these 
equations when temperature and or pressure 
effects are not significant produces negligible dif- 
ferences from the results at reference conditions. 



6. Conclusions 

High accuracy fluid measurements can be 
attained when the important factors affecting the 
performances of both flowmeter calibrators and 
fluid meters are properly taken into account in the 
measurement processes of these devices. To do 
this— to first order temperature and pressure 
effects on both the calibrator and the meter as well 
as on the fluid— the equations derived above 
should be used as the basic system models for cali- 
brator and meter. Specific features of calibrators 
need to be analyzed and, where pertinent, appro- 
priate modifications to the above-derived equa- 
tions should be made. Calibrators and metering 
units need to be properly instrumented and oper- 
ated according to the required assumptions and 
techniques for handling the data. Where further 
increases in flow measurement accuracy are 
needed, either more accurate descriptions of 
currently considered factors should be made or 
additional factors not currently considered should 
be assessed and included as pertinent, or both. 
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7. Appendix A. Cross-Sectional Area 
Changes in the Cylinder of the 
Calibrator 

The fluid pressure and temperature inside the 
cylinder of the calibrator can enlarge or reduce the 
cross-sectional areas of this cylinder, see Fig. 1. As 
well, geometrical changes of the rod or tube con- 
nected to the piston can also contribute to changes 
in this cross-sectional area. Of course, knowledge 
of speciflc geometries and material properties are 
crucial to characterizing these changes. In the fol- 
lowing it is assumed that a generic cylinder and rod 
geometry exist as sketched in Fig. 1. 

First, thermal expansion will increase the dia- 
meters of the cylinder and the connecting tube or 
rod via the first order approximations 



and 



D=Do[l + ac(Tc-To\ 
d=do[l + ar(Tr-To], 



(A.1.1) 
(A.1.2) 



where, in compatible units, D and d and £>o and do 
are the cylinder and tube diameters at non- 
reference and reference conditions, respectively. 
The linear expansion coefficients for cylinder and 
tube or rod are ac and a,, respectively, and Tc and 
Tt are the corresponding temperatures. These 
diametral enlargements will produce an increased 
cross-sectional area via the first order approxima- 
tion 



Ac = Aco [1 + 2ac (Tc- To] , 



(A.1.3) 



where Ac and Aco are, respectively, the annular ar- 
eas between cylinder and rod as averaged over the 
appropriate piston displacement. The linear expan- 
sion coefflcient of the rod or tube is assumed to be 
the same as that for the cylinder. 

Second, pressure effects can change the cross- 
sectional area through the enlargement of the cali- 
brator cylinder and the contraction of the rod or 
tube connected to the piston. The azimuthal or 
"hoop" stress, ag, produced in the calibrator cylin- 
der can be shown to be 



(T6 = (Pc-Po)Dc/2tc, 



(A.1.4) 



where Pc—Po is the pressure difference between 
the pressure in the cylinder and the reference pres- 
sure outside the cylinder. The quantities Dc and tc 
are, respectively, the inside diameter and the thick- 
ness of the cylinder. The azimuthal strain produced 
by this stress is assumed to be 



where transverse effects are neglected and where 
Ec is the modulus of elasticity of the material of 
the calibrator. This azimuthal strain increases the 
circumference of the calibrator cylinder via 



C=Co(l + ee), 



(A.1.6) 



where C and Co are, respectively, the inside 
circumferences of the cylinder in stressed and 
unstressed conditions. This produces an increase, 
to first order, in the cylinder's averaged cross- 
sectional area of 



' Ait Ait L 



{Pc-Po) gco l 
tcoEc J 



(A.1.7) 



or since Aa)= Co l^ir , 



^c=^co 1 + 



(Pc-Po) Dco l 



(A.1.8) 



Ce = aJEc , 



(A.1.5) 



where Z7co and ^co are the diameter and thickness, 
respectively, of the calibrator cylinder at reference 
conditions. Typical values for the ratio Deo/ too 
could be 10 to 20. 

In similar fashion, the cylinder pressure in excess 
of the reference pressure, Po at which condition the 
reference dimensions are determined can reduce 
the outside diameter of the tube connected to the 
piston. Since the deformation of a tube should, for 
the same material and diameter, exceed that for a 
solid rod, it shall be considered here. Therefore, 
the cross-sectional area of the tube connected to 
the piston would be 

- - r (Pc-Po) Dto^ 

where Dto and ^o are, respectively, the outside di- 
ameter and the thickness of the wall of the tube 
connected to the piston at reference conditions. 
The quantity £t is the modulus of elasticity of the 
tube material, at reference conditions. 

It is generally and safely assumed that the tube 
contraction effect expressed in Eq. (A.1.9) is negli- 
gible in comparison to the effects of temperature 
on the inner diameter of the calibrator. Therefore, 
it is assumed that the annular cross-sectional area 
between cylinder and tube is that expressed in Eq. 
(A.1.3) and the effect of pressure on this annular 
area is expressed in Eq. (A.1.8). However, where 
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speciflcs indicate that other compensations need to 
be made, other approximations can and should be 
done. The combined effect of pressure and temper- 
ature on the averaged cross-sectional area of the 
cylinder of the calibrator is, 



Ac=Acx>[l + 2ac 



(rc-rco][ 



n- 



iPc-Po)Dco 



Ectco 



]■ 



(A.1.10) 



8. Appendix B. Turbine Meter Data 
Processing 

Conventional practice in processing and using 
turbine-type flowmeter data began with plotting 
the meter K factor in pulses per volume units as a 
function of meter frequency, /; see Fig. A.2.1 and 
[5-8]. In this figure, the bracketed dots are 
intended to represent the mean values and the 
scatter obtained for the K factors at respective 
frequencies. For fluid and flow conditions which 
precisely duplicate those for which the meter cali- 
bration was performed, this practice produced 
acceptable results. However, when fluid or flow 
conditions deviated from those of the calibration, it 
was found that results could be unacceptable. An 
example is indicated in Fig. A.2.1. The lines for 
which the liquid viscosities are jm and jli2, respec- 
tively, indicate that different meter performance 
occurs at low flowrates when the fluid dynamic 
viscosity, /t2 is greater than the value jm used for 
the calibration. 



K 




Fig. A.2.1 Conventional format for processing turbine-type 
flowmeter data. 



Flowmeter performance which deviates from 
that of the calibration conditions at low flowrates is 
explained by citing the deviant phenomena which 



have, because of the different fluid or flow condi- 
tions, become signiHcant in comparison to the 
phenomena which prevailed during calibration. 
Shafer and Lee have resolved several situations 
involving viscous effects in specific flowmeter 
geometries [6,8]. From such efforts have come 
Universal Viscosity Curves (UVC) i.e., plots of 
meter factor versus the ratio of meter frequency to 
fluid kinematic viscosity. A typical UVC is shown 
in Fig. A.2.2. The ordinate, Ko, is the meter factor 
referenced to a selected reference temperature via, 
see [10], 



Ka=K[l+3aMiT-To]. 



(A.2.1) 



K, 




f/V 

Fig. A.2.2 Universal Viscosity Curve (UVC) for processing 
turbine-type flowmeter data. 



Deviations of flowmeter performance character- 
istics from UVC's can be explained by citing 
further deviants, such as lubricity, cavitation, 
extreme temperature or pressure effects on fluid 
properties, etc., see [5-8]. 

When the turbine meter is to be used in temper- 
ature conditions different from those of the refer- 
ence or calibration conditions, temperature effects 
may be compensated for by producing corrections 
for the meter factor K, as found in some standards, 
see [10]. This ISA standard produces an underived 
correction that is based upon the thermal expan- 
sion for the material of the meter body. This 
correction is effective in predicting trends for the 
case where the meter body and the turbine wheel 
material are the same. Where these materials are 
different, this correction can give erroneous 
results, for example, in the extreme situation where 
the turbine wheel has expanded to touch the inner 
wall of the meter body, thus stopping the wheel's 
rotation. In such a situation the meter factor 
should be zero yet the ISA correction would not 
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give this result. Therefore, further evolution is 
needed for conventional turbine-type flowmeter 
practice— beyond the UVC methods. 

Turbine meter performance is most properly 
characterized by specifying and quantifying, over 
the pertinent ranges of fluid and flow conditions, 
the significant dimensionless parameters which 
influence meter performance. To do this, the 
dimensional quantities that are involved in turbine- 
type flowmetering are cited as shown in Fig. A.2.3. 
These five dimensional quantities— fluid density p, 
and viscosity /x, the flow velocity U, the meter 
frequency /, and the meter diameter, D are 
described in terms of three independent units — 
length, time, and mass or force. Consequently, 
there are two dimensionless parameters to charac- 
terize the performance of this meter. These dimen- 
sionless parameters can take many forms. 
Conventional fluid mechanical practice would pro- 
duce parameters expressing ratios of meter 
frequency effects to fluid inertial effects, i.e., the 
Strouhal number and the fluid inertial to viscous 
effects, i.e., the Reynolds number. These are 
written 

St =fD/U (A.2.2) 

and 

DUp DV 
Re= — -= . (A.2.3) 




Fig. A.2 J Turbine meter configuration. 



However, turbine meters are designed to have 
frequencies proportional to volumetric flowrates. 
The constant of proportionality is the reciprocal of 
the meter factor. Since the constancy of this meter 
factor, i.e., the linearity of the meter over the 
flowrate (frequency) range is generally taken to be 
a measure of the quality of the meter, it has become 
conventional turbine meter practice to characterize 
meters via their K factor versus frequency depen- 
dence. In non-dimensional format, this can be done 
by specifying meter performance using a Strouhal 
number that is dependent upon the product of the 
Strouhal and Reynolds numbers, i.e., /D^/u. This 
parameterization is essentially that known as the 



UVC except for the factors involving the meter 
diameter, D, raised to an exponent. This product 
has been referred to as the Roshko number, Ro, 



Ro 



=!nl 



(A.2.4) 



see [9]. Alternatively, the pair of parameters, St 
and Re could also be used to characterize the 
meter's performance. In either case, only two 
parameters are needed to characterize the system 
sketched in Fig. A.2.3. 

In all of the above, the effect of different meter 
geometries or a different turbine wheel is not con- 
sidered. This is proper where the meter geometry 
is not changed and where the same turbine wheel 
and bearings, etc. are used and where the turbine 
wheel and the meter body are made of the same 
material. Under these circumstances, the above 
parameterization of the five dimensional quantities 
p,IM,U,D, and /is completed with the formulation 
of the Strouhal and Reynolds or Roshko numbers. 
Using this formulation the performance of the 
calibrated meter should enable the meter to be 
satisfactorily used over a range of conditions 
similar to those encountered in the calibration. 

If the situation is changed as shown in Fig. A.2.4 
so that the turbine wheel diameter, d, is included 
to produce six dimensional quantities, then three 
dimensionless parameters should be formed. The 
turbine wheels that are considered here are 
geometrically scaled versions of each other, i.e., 
they have the same number of blades, the same 
blade shape and only differ in the diameter, d. The 
nondimensionalization can be achieved by choosing 
the third dimensionless parameter to be /3=d/D, 
the ratio of turbine wheel to inside pipe diameter. 
The three parameter performance is sketched in 
Fig. A.2.5. It is noted that when the turbine wheel 
and meter body are the same material and when 
only thermal expansion effects are considered, p 
remains constant with changes in temperature, i.e.. 



P n 1 



D Doll + auiTM-ny 



(A.2.5) 



when Tt = TiA and at = au this diameter ratio is 
dolDo. 

When turbine wheel and meter body are differ- 
ent materials and when a,>au, the beta ratio is 
not constant and a condition for the linear expan- 
sion coefficients and the fluid flowrate can be 
derived for which Tu = T, and )3 = 1. This would 
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indicate that the turbine blade tips touch the inside 
of the meter body thus stopping the wheel. This 
condition is, of course, contrived in order to illus- 
trate the point that j8 variation can produce meter 
performance variation and in extreme limits can 
radically alter the graphs of meter performance. 




Fig. A.2.4 Turbine meter configuration with the turbine wheel 
diameter, d, considered as a parameter. 



St 



p 




Re or Ro 

Fig. A.2.S Non-dimensional turbine meter performance. 

If pressure effects are not negligible in altering 
the diameter of the meter body but are negligible 
in affecting the diameter of the turbine wheel, j8 is 
altered, to first order in temperature and pressure, 
separately, via 



St 



P(T) 



P=1 



Re or Ro 

Fig, A.2.6 Sketch of expected turbine meter performance where 

It is concluded that turbine-type flowmeter per- 
formance should be handled non-dimensionally. 
Pertinent corrections are then most clearly seen 
and most easily performed using the appropriate, 
first-order, relationships for the thermal and pres- 
sure effects. In this way high levels of flowmeter 
performance can be expected and achieved. 

9. Appendix C. Equations Summary 

The equations in Appendix C are numbered as in 
the main text. 

1. Calibrator Characterization : 

a. Geometrical Measurement Method— Refer- 
ence Conditions (see p. 513) 



Af Eo -Leo ^eo ^e 



A.CO — TT" — "=■ 



yea AcoLeu a 



(2) 



CO 



b. Geometrical Measurement Method— Non- 
Reference Conditions (see p. 515) 



P = 



do[i+ar(rr-ro)] 



Do [1 + au (Tm- To] 1 + — =- — ■=. 

L ^'MOiiM J 



(A.2.6) 

By so incorporating thermal and pressure effects 
using St, Re or Ro, and ^ parameters into perfor- 
mance curves as sketched in Fig. A.2.6, widely vary- 
ing metering conditions should be successfully 
handled. 



Kca—Kc 



where 



[l + aE(rE-ro)] 



[l-2ac(rc-r„][l ^^^^^J 






.(10) 



(7) 



c. Draw Technique— Reference Conditions 
(see p. 516) 



^CO— N^Eo/^'cOU.'TO . 



(15) 
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d. Draw Technique— Non-Reference Conditions 
(see p. 517) 



KcB-Kc 



{Pc-Po)Dco-] 



r (fc-fo)L)col 
[l + 2ac(rc-ro)][l+— ^^^^— J 



[1 -3aF(rcoLL'T- 7'c)][l + ^^^2!£2__^ jfj _ ^^(Te-To)] 



,(29) 



where 



Kc = Nn/ f'cOLL'T 



(28) 



2. Meter Calibration : 

a. Reference Conditions: (see p. 519) 



Data Processing: (see p. 522) 
St =fuD ulUu^ CiKuD M , 

Re=DuUulv , 

„ fuDu 
Ro=—^ • 

3. Flow Measurement: 

a. Reference Conditions: (see p. 522) 



• /mo CiZ) mo/mo 
Kmo— zi — = := . 



b. Non-Reference Conditions: (see p. 523) 



(54) 
(55) 

(56) 



(59) 



/mo -Keo ^Vmo -Keo Nii 



/m 



^Mo-7 5=~-1\J =^--Tj— •-K'co=7-*-S^ca. 

/EO y4cO •'^EO ^co •'^EO /EO 

(42) 



Dflto Processing: (see p. 519) 



5/— — TT —CiKmqDmo^KmoDmo, (43) 



f/ 



MO 



Re = 



D mhU MO 



Ro 



/mo-Dmo 



(44) 



(45) 



b. Non-Reference Conditions: (see p. 521) 



'tM-7"Ac -T7~"-C0 

/e PC JVe 



[1 - gE(rE-7'o)] [1 + 3aF(7'c- Tm)] 



(53) 



• _/m _Ci/m-Dm 



^N 



5r 



(60) 



, Wl,.«o,f>--(--^-*^^^l^] 



5/ ^M^M 






[l+3ap(rM-7'Mo)][l-^^^] 



(65) 
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